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I. INTRODUCTION 

The properties of strongly-interacting matter at non- 
vanishing temperature and density are in the focus of 
many theoretical and experimental efforts. On the exper- 
imental side, future and running heavy-ion experiments 
at various facilities such as GSI, JINR, CERN and BNL 
aim at probing the phase structure of quantum chromo- 
dynamics (QCD), especially in the regime where poten- 
tially a critical endpoint (CEP) is present. 

On the theoretical side, functional continuum methods 
are well-suited for a combined study of the confinement- 
deconfinement and chiral phase structure of QCD at fi- 
nite density. In recent years, functional renormalisa- 
tion group (FRG) studies and Dyson-Schwinger studies 
have provided valuable insights into the phase structure 
of strongly interacting matter. These studies were per- 
formed both within first principle QCD as well as within 
low-energy effective models, see e.g. [1-7]. For QCD- 
related reviews see [4, 8-16] (FRG) and [17-22] (DSE). 

At vanishing density, lattice simulations provide addi- 
tional insights from first principle QCD that are com- 
plementary to the continuum studies. This also helps to 
improve the systematic error analysis of the respective 
approaches. At non-vanishing chemical potential, how- 
ever, lattice simulations are hampered by the sign prob- 
lem, see e.g. [23, 24]. Impressive progress in overcoming 
this problem has been made [25-31], but to date lattice 
QCD is still restricted to small chemical potential. 

Low-energy effective models, such as the Nambu-Jona- 
Lasinio (NJL) and quark-meson (QM) model, are known 
to describe the chiral dynamics of QCD rather well, for a 
recent review see e.g. [32]. Polyakov-loop extended ver- 
sions thereof (PNJL, PQM) [33, 34], have advanced our 
understanding of the confinement-deconfinement aspects 
of the QCD phase structure, see e.g. [2, 3, 33-44]. It has 
been known for a long time, see e.g. [9, 11, 12, 14, 15], 
that low-energy effective models can be systematically 
related to full QCD within the FRG approach. For 
Polyakov-loop extended models this follows from the 
Landau gauge approach in [1, 4, 45], for the Polyakov 
gauge see [46, 47]. This setting has been discussed in de- 



tail for the PQM model in [3, 4, 48], see also our discus- 
sion in Sec. II. While chiral symmetry and its dynamical 
breaking are well described within these models, confine- 
ment is only included in a statistical manner. Moreover, 
the glue potential of full QCD, encoding the gauge dy- 
namics in the presence of matter fields, is replaced by a 
Polyakov-loop potential. This potential is fixed to lattice 
data of the pure Yang-Mills system at vanishing chemical 
potential. In such an approach, the coupling of the mat- 
ter sector to the gauge sector is lost, see [3, 34]. This is 
discussed in detail below, where we show how the matter 
back-coupling can effectively be taken into account. 

In conclusion, the FRG-approach to the low-energy dy- 
namics of QCD is a promising set-up for investigating 
the QCD phase structure and thermodynamics. A study 
along these lines has already been put forward in [3, 48]. 
In the present work we extend this approach by a more 
refined inclusion of the matter and glue dynamics as well 
as studying the mass sensitivity of various observables. 

This work is structured as follows: In Sec. II we intro- 
duce the Polyakov-quark-meson model as a low-energy 
truncation of full two-flavour QCD. Furthermore, the in- 
clusion of the matter back-coupling to the gauge sector 
is discussed. Sec. II A summarises the renormalisation 
group approach and the resulting flow equation for our 
model. In the following Sec. Ill we compare the phase 
structure of the PQM model to that of the QM model 
and discuss the mass sensitivity of the phase transitions. 
In particular, we find an intriguing splitting of the chiral 
phase transition at low Goldstone-boson masses, which is 
only observed when fluctuations are taken into account. 
Some thermodynamic observables and their mass depen- 
dence are studied in Sec. IV. Concluding remarks and 
an outlook are given in Sec. V and some details of our 
numerical implementation can be found in the Appendix. 



II. FROM QCD TO POLYAKOV-LOOP 
EXTENDED CHIRAL MODELS 

We have already emphasised in the introduction that 
low-energy effective models can be related to full QCD in 
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FIG. 1: Partially hadronised version of the FRG flow for 
QCD. The loops denote the gluon, ghost, quark and meson 
contributions, respectively. The crosses mark the RG regula- 
tor term. 



a systematic fashion. In [1, 3, 4, 46-48] this is discussed 
for the embedding of Polyakov loop extended models: 
The FRG provides a setting which allows to approach the 
low temperature regime of QCD from large momentum 
scales by successively integrating-out momentum shells. 
This entails that one starts at asymptotically large mo- 
menta with pcrturbative QCD with one coupling param- 
eter, the strong coupling as(p^). The other parameters 
are the entries in the quark mass-matrix. In the present 
two-flavour study the masses of up and down quark are 
small and the mass matrix is diagonal in flavour space. 

When lowering the momentum scale one systematically 
includes quark and gluon fluctuations into the theory, fi- 
nally approaching the hadronic phase, see [49] for one 
flavour QCD and [1, 4] for two-flavour QCD. Close to the 
phase boundary between the quark-gluon plasma phase 
and the hadronic phase and at not too large chemical po- 
tential, mesonic degrees of freedom, in particular the pion 
and sigma fluctuations, become important. The related 
quark bilinears, qq and qjsfq, carry the same quantum 
numbers, i.e. they have a considerable overlap with the 
full meson operators. The matter sector of the QCD ef- 
fective action can be conveniently written in powers of 
these bilinears and related kinetic terms. The strength 
and the momentum dependence is systematically com- 
puted by means of the flow equation for QCD, depicted 
in Fig. 1. The propagators in Fig. 1 are the fully dressed 
field dependent propagators and the flow equation de- 
scribes the fully coupled QCD glue-matter system. In 
particular, there are contributions of the matter sector 
also in the diagrams for the gauge sector, i.e., the vac- 
uum polarisation diagram in the gluon propagator, see 
Fig. 2. The first two loops in Fig. 1 constitute the glue 
potential of QCD. Note that this is not the Yang-Mills 
potential, as the gluon and ghost propagator are those of 
QCD and scale differently with momenta. In turn, drop- 
ping the matter back-coupling these two loops correspond 
to the flow of the glue potential of pure Yang-Mills the- 
ory [50-52], for more details see [4, 53]. The remaining 
two loops represent the flow equation of the dynamical 
quark-meson model, see e.g. [54, 55]. 

Following these arguments, we employ the Polyakov- 
quark-meson (PQM) model [3, 34] as an effective realisa- 
tion for low-energy QCD. As mentioned in the introduc- 
tion, the PQM model is an extension of the quark-meson 
(QM) model by Polyakov-loop variables which allows to 
effectively describe chiral and certain aspects of confine- 
ment in QCD. The Euclidean Lagrangian of this model 



FIG. 2; Quark polarisation contribution to the gluon propa- 
gator representing a contribution to the matter back-coupling. 



for Nf — 2 flavours and Nc — 3 colours including a uni- 
form quark chemical potential /i, reads 

-CpQM = + /i(o- + i75T7f) + A*7o) g + ^(<9pcr)^ 

+ \id^7Tf + U{a,7r)+U{<f,^;To), (1) 

where (j) = (cr, tt) denotes the 0(4)-symmetric representa- 
tion of the meson fields. The scalar sigma meson and the 
three pseudoscalar pions are coupled via a flavour-blind 
Yukawa coupling h to the quark fields, q and q. The co- 
variant derivative = 7^9^ — i g7oAo(*&) couples the 
Polyakov loop, defined as 



^^tr7'exp|^ig^^dTAo(f,T)^^ , (2) 



to the quark fields. In this work we describe the theory 
solely in terms of the Polyakov loop $ and suppress the 
Ao-dependence. Differences between these two formula- 
tions will be discussed elsewhere. Moreover, we add a 
purely mesonic potential ?7(cr, 7?) that allows for explicit 
as well as spontaneous breaking of chiral symmetry 



A, 



(3) 



The explicit breaking of chiral symmetry is realised by 
a term linear in the sigma field. The pure quark-meson 
(QM) model is obtained from the PQM model by setting 
Aq = 0, i.e. $, $ = 1 and omitting the glue potential 
W($,$;ro). 

In Polyakov-loop extended chiral effective models, the 
glue potential of full QCD is usually approximated by 
an effective Polyakov-loop potential i>; Tq) which is 
constructed by means of a Landau-Ginzburg-like ansatz. 
The arising coefficients are fitted to lattice results for the 
pure-glue system, see e.g. Eq. (4) below. The explicitly 
shown parameter Tq determines the scale of the decon- 
finement phase transition. 

In the following we apply a simple polynomial ansatz 
for the glue potential as introduced in [39] 



rp4 



with temperature-dependent coefficient 



(5) 
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The parameters of Eqs. (4) and (5) have been determined 
in [39] by a fit to pure Yang-Mills lattice results to be 



flo = 6.75 , ai 
and 



-1.95 , 02 = 2.625 , aa = -7.44 (6) 



&, = 0.75 , &4 = 7.5 . 



(7) 



In this construction, however, the back-reaction of the 
matter sector is completely neglected. This deficiency 
can be overcome by replacing the constant parameter Tq 
by a flavour number, Nf, and quark chemical potential, 
/X, dependent function: Tq To{Nf,T,fi). For the Nf- 
and //-dependence of Tq we follow the arguments in [3, 
34], to wit 



with the T-scale and 



(8) 



6(iV/,T,/i) 



- b„, ,„ Bt,{T,^j.) , (9) 



67r 



with niq — h (Tvac and CTvac = /tt = 93 MeV. For details 
on our choice of parameters see [3]. The /i-dependent 
correction in Eq. (9) is proportional to a difference of 
(baryonic) Fermi-Dirac distributions 

0,(r, ^i) = UhiT, fi) + n,(r, -/i) - 2nb(r, o) , (lo) 



with 



nb(T,fi) 



1 

1 + e3im^-f^)/T 



(11) 



The physical importance of this expression is seen in the 
limit of vanishing temperature. Then 9i, reduces to the 
step function — m^), as it should in order to account 
for the Silver-Blaze property of QCD: at vanishing tem- 
perature, chemical potential effects are expected to only 
set in when /i > rUg. Note that in the present approxi- 
mation the binding energy of the nucleons is neglected. 
As the temperature increases, the sharp behaviour of the 
step function is smoothed out. 

The origin of b{Nf , T, fj,) is the flow of the fermionic 
part of the vacuum polarisation, see Fig. 2, or more pre- 
cisely 



(12) 



in the presence of the Polyakov loop. This term car- 
ries the fermionic part of the QCD /3-function. Its 
//-dependence stems from the Polyakov-loop enhanced 
quark/anti-quark occupation numbers 

7V,(r,/i;$,$) (13) 

~ 1 + 3<|e(£^g-A')/r ^ 3$g2(£;,-M)/T g3(B,-At)/T ' 



(fc2 



igV^^, which has been used in [1, 4]. 



with Eg 

Note that the occupation number Eq. (13) appears in 



any quark loop with a Polyakov loop or Aq background. 
In fact, the flow of the fermionic part of the free energy 
both in QCD, [1, 4] and in the PQM-model, [2, 3] as well 
as Eq. (15), are proportional to Eq. (13). The latter flow 
is also central for the present work, see Eq. (15) below. 

Here we use Eq. (13) to estimate the //-dependence 
of the QCD /3-function at vanishing cutoff scale, fc = 0: 
For large temperatures T > Tc{^) and hence deconfining 
Polyakov loops, $,$—)■ 1, the thermal distribution factor 
in Eq. (13) reduces to the standard Fermi-Dirac distribu- 
tion for quarks, see also [49]. In this large temperature 
limit we also infer that Ng — >■ 1/2. In turn, for tem- 
peratures T < Tc{fJ-) and confining Polyakov loops with 
$, $ — > 0, and a = (Tvac, we arrive at nfc(T, /() in Eq. (11). 
Due to the $, ^-independence of the large temperature 
limit we simply use confining Polyakov loops $, $ = 
for all temperatures. This leads to our final expression 
Eq. (10). 

In summary, the non-trivial factor Of, in the definition 
To{Nf,T,fi) improves the original choice of 6(7V/,/i) in 
[34] with respect to the thermodynamical properties by 
taking into account the dynamical change of the quark 
contributions to the glue sector for finite temperature 
and chemical potential. 



A. Flow Equation for the PQM Model 

For a realistic description of phase transitions one has 
to include thermal and quantum fluctuations. In the 
present work this is done within the functional renor- 
malisation group (Wetterich) equation [56] 



dtTkix] = ^STr 



Ri 



dtRk 



(14) 



In this equation t — log(/e/A) denotes the RG-time, x 
represents a generic field variable and Rk is the RG reg- 
ulator that implements the idea of integration over mo- 
mentum shells. The super-trace STr, involves a trace 
over internal (colour, flavour and Dirac) spaces as well as 
a momentum integration. Furthermore, it accounts for 
the correct signs of the fermionic and bosonic contribu- 
tions. 

After integrating-out the gluonic degrees of freedom, 
the flow of the free energy contains solely the last two 
loops in Fig. 1, while the first two lead to the glue po- 
tential rigiuc($, 5')- Moreover, this procedure results in 
modifications of the matter sector due to the coupling to 
the Polyakov loops. This is accounted for by our initial 
conditions for the quark-meson sector. The QCD free 
energy is thus given by 

f^QCD(CT,7r,$,l>) = r2giuo($,$) +anattor,A((T,7r,$,$) 



dfc(9fcl]i„attor,fe(CT,7r,$,$). (15) 



Since the gluons have been integrated out, their dynamics 
is stored in the Polyakov-loop glue potential ficriuo- More- 



4 



over, the full QCD free energy is necessarily independent 
of the cutoff, i.e., i9a^qcd = 0. Hence, A-dependent 
terms appear in rimattor.A- These can be determined 
from the flow at A, see e.g. [8, 9, 11-13, 15]. From a 
phenomenological point of view these contributions cor- 
respond to the high-energy part of the vacuum fluctua- 
tions [57]. Finally, the flow equation for the matter sector 
of two- flavour QCD in a Polyakov loop or ^o-background 
reads [1-3] 

d,n, = ^ / ^ coth (^] + ^ coth ( ^] 



;i-7V,(r,/.;$,$)-iV,(r,-/i;$,$) 



(16) 



with the Polyakov loop enhanced particle numbers Nq 
defined in Eq. (13). The quasi-particle energies are given 
by Ei — y/k!^ + mf, i = q, 7r,(T, and the masses are de- 
fined as 



m 



2 J,2 



= 2n[ 



k I 



(17) 



In the above expressions, a prime at the potential denotes 
the derivative with respect to 0^. 

The order parameters xo = {<^0i ^o) for given tem- 
perature and chemical potential arc determined by the 
solution of the corresponding equation of motion (EoM) 



da 



fe-i-O 



9$ 



5$ 



. 



(18) 



For the numerical solution of the coupled Eqs. (18), 
we utilise a stochastic technique which is outlined in 
App. A. In comparison to standard multi-dimensional 
root-finding algorithms, such as Newton's method, we 
obtained with this technique a much higher numerical 
accuracy within reasonable CPU time, which is required 
in particular for the evaluation of thermodynamic quan- 
tities. 



B. Initial Condition at Vanishing and Finite 
Temperatures and Density 

In order to solve the flow equation (16) numerically, the 
parameters A, v'^ and c of the meson potential U{<7, n) in 
Eq. (3) as well as the Yukawa coupling h have to be speci- 
fied at the UV scale. However, the parameters are not in- 
dependent and are related to vacuum low-energy observ- 
ables in the IR. For example, the explicit chiral symme- 
try breaking parameter c relates the pion decay constant 
and the pion mass via m^/^ = c and can thus be fixed 
by these observables. The remaining three parameters 
at a given UV cutoff are chosen such that specific values 



of low-energy observables (/jr, Too-, rUq) are reproduced in 
the IR. In particular, the physical mass point is char- 
acterised by the values = 93 MeV, to^ — 540 MeV, 
rUq = 297 MeV and = 138 MeV. To achieve these 
values in the infrared we have chosen the following initial 
values at the UV scale A = 950 MeV: 



A 

c 
h 



1.3, 

-4.36 • 10^ (MeV)^ , 
1.77 • 10^ {MeYf , 
3.2. 



(19) 



Implicitly, Eq. (19) also takes into account the gluonic 
fiuctuations: below the UV scale A they effectively decou- 
ple, their contributions above the UV scale lead to the ini- 
tial conditions Eq. (19) and the Polyakov-loop potential. 
Moreover, for sufficiently large UV scale A, thermal and 
chemical potential modifications are negligible as they 
are suppressed exponentially with — A/T for thermal fluc- 
tuations and polynomially with /i/A for density-related 
fiuctuations. For large temperatures and/or chemical po- 
tential, however, the initial conditions in Eq. (19) receive 
corrections. 

In general, thermal and quantum fluctuations with mo- 
menta A: < A have to be taken into account in a tempera- 
ture and chemical potential dependence of Eq. (19). This 
can be either done by solving the related QCD flow for 
fc > A or by projecting this flow on the matter sector. 

In our previous work [3] we have improved the thermo- 
dynamic observables by including the integrated UV-flow 
in the initial condition of the free energy. 



f^mattor,A [cr, 7?, $] = C/(cr, 7?) 



17^[a,7f,$,$], (20) 



evaluated on the EoM. Here, fl"]^ corresponds to the in- 
tegration over the flow for the interacting Polyakov-loop 
system for scales fc > A and fixed parameters. This en- 
sures the A-independence of system at k = 0, i.e. RG- 
invariance. A more detailed discussion this procedure 
can be found e.g. in Refs. [2, 3, 57]. 

In the present work we extend this procedure to the full 
initial condition of the effective action. This is implicitly 
based on the results for the coupling flows in QCD with 
partial or full dynamical hadronisation [1, 4, 58]: the flow 
of the QCD coupling parameters directly related to the 
parameters Eq. (19) is small, and hence the computation 
of the thermal and chemical potential modifications for 
fc > A on the basis of the T — values at A is already 
a good estimate of the full modification in QCD. More 
details on this will be presented elsewhere. 



III. MASS SENSITIVITY OF THE PHASE 
STRUCTURE 

For the subsequent analysis of the mass sensitivity of 
the phase structure and thermodynamics, we will only 
vary the explicit chiral symmetry breaking parameter c. 
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which determines the Goldstone-boson mass and keep all 
other parameters fixed. 

From the order parameters we can deduce the chiral 
and deconfinement phase transition lines T^{fi), Td{^), 
respectively, in the (T, /i)-phase diagram. At physical 
pion masses and small chemical potential, both transi- 
tions are crossovers, entailing that there exists no unique 
definition of the transition temperature. In the following 
we use the inflection point of the corresponding order 
parameter to define a (pseudo-) critical temperature. In 
addition, we compare this transition point to the one de- 
fined by one half of the normalised order parameter 

a(0,0) "2- ^^^> 

In order to highlight not only the influence of fluc- 
tuations, but also the direct impact of the Polyakov 
loop on the phase structure, we compare results obtained 
with the PQM model with those of the pure QM model. 
Furthermore, all PQM results presented below include 
the matter back-coupling To(7V/, T, /x), as introduced in 
Eq. (8). 



A. Physical Mass Point 

We begin with a discussion of the phase structure for 
physical values of the low-energy observables. The cor- 
responding (T, /i)-phase diagram for the QM model is 
shown in the left panel of Fig. 3. At low quark chemical 
potential we find a chiral crossover around T k, 160 MeV 
which turns into a first-order phase transition at large 
chemical potential in a critical endpoint (CEP). 

In comparison to standard mean-field results, where 
meson fluctuations are ignored, we observe that the lo- 
cation of the CEP is shifted towards lower temperatures 
when fluctuations are taken into account, see e.g. also, 
[55]. Interestingly, a similar behaviour is found if the 
standard mean-field approximation is improved by con- 
sidering the renormalised QM models. The renormalisa- 
tion of these models amounts to including vacuum fluc- 
tuations of the quark loop [59, 60]. Already these fluc- 
tuations push the location of the CEP towards higher 
chemical potential and smaller temperatures and exclude 
the existence of a CEP at low fi/T ratios. A similar con- 
clusion can be drawn for renormalised (P)QM models 
with three quark flavours [61]. Note that these vacuum 
fluctuations are always included in the FRG treatment. 

In addition, in Fig. 3 we compare the transition line 
obtained by the inflection point of the chiral order pa- 
rameter (short-dashed line) with the one obtained via 
Eq. (21) (long-dashed). At low chemical potential both 
curves agree well, but start to deviate at larger chemi- 
cal potential . The curve defined by the inflection point 
runs into the CEP and the transition turns into a first- 
order one for smaller temperatures. A for RG calcula- 
tions typical back-bending of the transition curve at low 



temperatures towards smaller chemical potential is ob- 
served. This behaviour is not seen for the other defini- 
tion of the chiral transition line, which does not reach 
the endpoint. In contrast, the curve bends outwards at 
lower temperatures and hits the /i-axis at /i « 322 MeV. 

In the right panel of Fig. 3 we include the Polyakov 
loop and show the PQM phase structure for the same 
IR values. At vanishing and low chemical potential, the 
chiral and deconfinement transitions are both crossovers 
and again both definitions of the chiral transition are 
shown. For the chiral transition, a similar behaviour as 
in the pure QM model is found. 

The Polyakov-loop related transitions, which we refer 
to as "deconfinement" transition lines, are denoted by 
the two black (dot-dashed and dot-dot-dashed) lines, de- 
fined by the inflection points of the Polyakov loop and 
its conjugate. The dark band in Fig. 3 denotes the width 
of d^/dT at 80% of its maximum height and measures 
the strength of the phase transition. At vanishing chem- 
ical potential the chiral and deconfinement transitions 
lie close to each other and continue to coincide for in- 
creasing chemical potential. However, if the matter back- 
coupling is neglected, a splitting of the two transitions at 
finite chemical potential has been observed in previous 
mean-field [34, 61] and RG calculations [3]. In partic- 
ular, the deconfinement transition line becomes almost 
/i-independent , i.e., ^^(/i) ~ Td{^ ~ 0) without the mat- 
ter back-reaction. This effect supported previous specu- 
lations about a possible quarkyonic matter region in the 
QCD phase diagram, where chiral symmetry is restored 
while deconfinement persists [62]. In the present two- 
flavour study including fluctuations as well as the matter 
back-reaction we, however, do not observe such a region 
anymore. A similar scenario has also been found in re- 
cent non-perturbative Dyson-Schwinger studies for two- 
and (2 -f l)-flavours, cf. [7, 63] and references therein. 

Towards the critical endpoint, a narrowing of the dark 
band is seen that indicates a sharpening of the transi- 
tions. For temperatures below the CEP, the chiral tran- 
sition is of first order. Also in the PQM model, the 
CEP is located at low temperatures, (a^cePj Tcep) ~ 
(290, 20) MeV. 

At this point, a word of caution concerning the high- 
chemical potential region should be added. In this re- 
gion important baryonic degrees of freedom and diquark 
fluctuations are usually ignored in the literature. Due 
to their complexity this is also done in the present work. 
Recent estimates of the influence of baryons on the phase 
structure for two colours can be found, e.g., in [64-67]. 
Due to these omissions the phase structure at high chem- 
ical potential becomes questionable. However, the region 
that is to date accessible to lattice calculations, [ijT < 1, 
is well within the region of applicability of the present 
study. Hence, we can safely conclude that for low chem- 
ical potential the existence of a critical endpoint can be 
excluded, see also [68]. 
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FIG. 3: Phase structure at physical pion masses for the QM (left) and PQM (right) models. Depicted is the chiral (blue, long- 
and short-dashed) transition line. In the PQM model the Polyakov-loop transitions (black, dot-dashed and dot-dot-dashed 
lines) are also shown. 



B. Chiral Limit 

In the following, we investigate the mass sensitivity of 
the phase structure. As an extreme case we consider the 
chiral limit, where no explicit chiral symmetry breaking 
terms are present and the Goldstone bosons, the pions, 
are massless. In this limit the chiral transition at van- 
ishing chemical potential is a sharp transition of second 
order lying in the 0(4)-universality class, cf. [69]. The 
resulting chiral phase structure is shown in Fig. 4. Com- 
pared to physical masses the critical temperature of the 
chiral transition decreases by approximately 20 MeV due 
to the decreasing constituent quark masses. Since there 
is no ambiguity in the definition of the transition tem- 
perature in this limit, we only show one transition line. 

At high chemical potential and low temperatures, a 
novel phase structure emerges: the chiral transition line 
splits into two branches. This behaviour has previously 
been observed in a QM model RG study [54], see left 
panel in Fig. 4 for our QM result. In accordance with 
the previous findings we find two critical points (CPs) in 
this region: one on the outer transition branch at low T 
and a further CP at higher temperature. 

However, for the PQM model only one CP survives 
which is located on the inner branch of the transition 
line, at (/icPjTcp) = (255,27) MeV, similar to the loca- 
tion for physical masses. Below this point, the transition 
is of first order, while the outer branch remains of sec- 
ond order. Note that due to the finite quark masses, 
the Polyakov-loop related transitions are still crossovers. 
Furthermore, the sharp transition in the chiral sector in- 
duces an additional peak in the temperature derivatives 
of the Polyakov loops. This interferes with the definition 
of the dark band in the PQM phase diagram and sub- 
sequently, the focusing of this band towards the critical 
point is not directly seen in this limit. 

Furthermore, the chiral splitting is also reflected in the 
behaviour of the sigma meson mass as a function of the 



chemical potential as illustrated in Fig. 5 for the PQM 
model at T = 30 MeV, i.e. just above the critical point. 
At low temperatures two minima, corresponding to the 
two branches of the chiral transition, appear in m(j(/i). 
In contrast to the sigma mass, the pion mass remains 
zero until chiral symmetry is completely restored, which 
happens at the outer transition branch. The quark mass, 
on the other hand, is proportional to the chiral order pa- 
rameter, niq = ha. Thus it shows a change in slope at the 
first transition and reaches zero at the second transition. 



C. Small Pion Mass 

Here we demonstrate how the splitting region in the 
phase diagram changes when the pion mass is increased 
towards the physical point. 

The resulting phase structure for — 50 MeV is 
shown in Fig. 6. Due to the non-vanishing explicit chi- 
ral symmetry breaking, the chiral transition immediately 
becomes a crossover. In the low temperature/high chem- 
ical potential region we still observe a splitting in the 
chiral transition line and also the second minimum in 
the sigma meson mass persists. As the pion mass is fur- 
ther increased towards its physical value, however, the 
outer branch of the transition is weakened. 

We also comment on the impact of the thermal and 
chemical potential modification of the initial conditions 
discussed in Sec. II B on the above effects: When the ini- 
tial conditions of the flow are fixed in the vacuum and 
T-independent To(/i) is used, it is not possible to define 
a second branch of the chiral transition at the physi- 
cal point. The inclusion of a T and fj, dependence of 
the initial action, Ua{T,ij,), as well as the more sophisti- 
cated definition of To(T, fj,) in Eq. (8), however, enhances 
this effect once more, as shown in Fig. 7. In this figure 
we compare the /x-derivative of the chiral order param- 
eter for To(/i) and rJA(T,^) = 1^a(0,0) ("standard") to 
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limit as a function of /i for constant temperature T = 30 MeV. 
One extra minimum of rria occur in the splitting region while 
vanishes until the second chiral transition is reached. 



the result using To(T, /i) and Vlf^iT^fji) ("enhanced") at 
physical masses. While both versions show a peak at 
/i^ « 293 MeV - defining the (inner) transition branch - 
the enhanced version shows a clear additional peak cor- 
responding to the second branch. In fact, the line defined 
by this second peak agrees well with the one determined 
by Eq. (21) and is hence not explicitly shown in Fig. 3. 

Note also that the location of the critical point in tem- 
perature direction, Tcep, is only mildly sensitive to the 
pion mass. In all considered cases we find a critical 
point at TcEP ~ 20-30 MeV. This is in contrast to pre- 
vious PQM model results where (Tcep, Mcep) depends 
strongly on rria when ttItt , f-K and ruq were fixed to their 
physical values, see [70] for mean-field results. When fiuc- 
tuations beyond the mean-field approximation are taken 
into account, this observation still holds. Here, however, 
we vary only the explicit chiral symmetry parameter, 
which results in different values for e.g. and m^. 

On the other hand, the location of the critical point 
along the chemical potential axis, /icEP, is changed more 



drastically as is varied. This can be understood by 
noting that the critical chemical potential at vanishing 
temperature, iJic{T = 0), is related to the quark Fermi 
surface and hence to the quark mass which increases with 
m^r, and so does ^ic{T = 0). 



IV. THERMODYNAMICS 

In order to achieve deeper insights into the nature of 
the phase transitions and the impact of the Polyakov 
loop, we focus in the following on some thermodynamic 
observables. 



A. Thermodynamic Observables 

The thermodynamic grand potential is obtained by 
evaluating the infrared effective average potential on the 
EoM (18) 



n(T,^l) = fife^o(r,M)l 



(22) 



We begin with a discussion of the pressure which is de- 
fined as the negative value of the grand potential 



p{T,^x) = -^{T,^i) + ^{o,{)) 



(23) 



and is normalised to zero in the vacuum. Since the pres- 
sure represents a thermodynamic potential, all further 
observables follow from this expression in the standard 
way by differentiation. For example, the first derivatives 
of the pressure with respect to T and ^ yield the entropy 
and quark number density 



gp(T,/z) 
dT 



(24) 



respectively. For high temperatures and densities we use 
the Stefan-Boltzmann pressure of massless QCD, corre- 
sponding to an ideal gas of massless quarks and gluons 
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(25) 



The expression in the first line denotes the fermionic con- 
tribution and the second line contains the gluonic part. 
Appropriate derivatives of this expression can be used to 
normalise other thermodynamic observables. 

In addition wc can define the energy density e = —p + 
Ts + 2/irig which is used to calculate the trace anomaly 

A _ e;; _ e - 3p 

rp4 



ji4 



rp4 



(26) 



which is related to the trace of the energy-momentum 
tensor 0'^" and vanishes in a scale invariant theory. This 



quantity thus yields a measure for the breaking of confor- 
mal invariance of the system. Furthermore, this observ- 
able is also referred to as the interaction measure, since 
it quantifies the deviation from the equation of state of 
an ideal gas e = 3p, cf. Eq. (26). 



B. Mass Sensitivity 

In the remaining section we investigate the thermody- 
namic observables, the impact of the Polyakov loop on 
these in hot and dense matter and study their mass sen- 
sitivity by varying the pion masses. 

Fig. 8 shows the pressure normalised by its Stefan- 
Boltzmann value for three fixed values of chemical po- 
tential and pion masses ttZtt = 138, 50, MeV from left 
to right. The largest chemical potential is chosen such 
that we pass close to the critical point. Results for the 
QM model are depicted in the upper panels, while the 
PQM calculation is presented in the bottom panels. 

The QM pressure levels at p/psB ~ 0.6, owing to the 
lack of gluonic degrees of freedom. At low tempera- 
tures, the pressure is dominated by the lightest mesons, 
i.e. pions. Due to the complete lack of confinement in 
the QM truncation, the pressure rises almost immedi- 
ately, while the slope becomes steeper for higher chemi- 
cal potential. As the pion mass is lowered, the expected 
dominance of the Goldstone bosons in p is visible. In 
the chiral limit and at vanishing chemical potential, a 
plateau develops that is slightly below the expected value 
of p/T^ = 37r^/90 for a free gas of massless pions [71], 
denoted by the black line at low T. 

Including the Polyakov loop, which represents a statis- 
tical implementation of confinement, the picture changes. 
Now, the pressure is almost constant in the hadronic 
phase. Again, a plateau develops at low temperatures 
that becomes more and more pronounced as we decrease 
m^. With increasing temperature, as the quark masses 
decrease due to the restoration of chiral symmetry, the 
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FIG. 8: Pressure normalised by the Stefan-Boltzmann pressure, Eq. (25), in the QM (upper panels) and PQM (lower panels) 
model, each for — 138,50,0 MeV (left to right). 



pressure rises more strongly. As the chemical potential 
increases, the chiral and deconfinement transition tem- 
peratures are lowered and the pressure rises strongly al- 
ready at lower temperatures. It saturates at about 80- 
90% of the Stefan-Boltzmann limit at high temperatures. 
A similar behaviour is found for all pion masses. 

Decreasing the pion mass, however, another problem 
arises: The coupling of the Polyakov-loop and quark sec- 
tors seem to be badly balanced. As a result, we find that 
the pressure increases slower than at /.t = 0, resulting 
in a non-monotonic ratio p/psB- We emphasis that nev- 
ertheless, p itself is a monotonically increasing function 
of temperature, as it should be. 

Combining the energy density and the pressure as de- 
fined in Eq. (26), we obtain the interaction measure A/T^ 
which we show in Fig. 9. This quantity again illustrates 
nicely the improvement achieved by the inclusion of the 
Polyakov loop. The QM result of this observable showed 
an unphysical two-peak structure at non- vanishing chem- 
ical potential. This effect is cured by the inclusion of 
gluonic degrees of freedom. The interaction measure is 
approximately zero at low temperatures and increases 
strongly around the phase transition. At high tempera- 
tures it decreases ~ l/T^. Of course, the non-monotonic 
behaviour of p/T^ at < 138 MeV also influences this 
quantity, which results in slightly negative scale anomaly 
at low T for these masses. For physical masses, however, 
we even find good agreement with the recent lattice re- 
sults, cf. [43]. 



V. CONCLUSIONS AND OUTLOOK 

Based on a renormalisation group analysis we have 
presented results on the phase structure and thermody- 
namics of the Polyakov-quark-meson truncation for two 
quark flavours, which serves as an effective model for low- 
energy two-flavour QCD. Special emphasis is put on the 
influence of thermal and quantum fluctuations, that are 
crucial for a proper description of phase transitions. In 
particular, we have argued that the PQM model consti- 
tutes a well-controlled approximation to first-principles 
full QCD. One focus is given to the back-coupling of 
quarks to the glue sector of QCD which results in a Nf-, 
T- and /x-dependent modification of the To parameter in 
the Polyakov-loop potential. 

At physical pion masses we find a chiral critical end- 
point at much lower temperatures than found in standard 
mean-field calculations. This nicely demonstrates the im- 
portance of fiuctuations. Moreover, we can rule out the 
existence of a critical point for small chemical potential 
^/r« 1-2. 

When studying the mass sensitivity of the chiral phase 
diagram, we encounter an intriguing phase structure. In 
the chiral limit, at high chemical potential and low tem- 
peratures a splitting of the chiral transition line into two 
branches is observed. This splitting is not observed in 
the mean-field approximation in the same model, even if 
the vacuum term of the quark loop is taken into account. 
A similar effect has, however, previously been found in a 
FRG study of the quark-meson model [54]. Differences 
between the previous and the present work can be at- 
tributed to our higher value of the UV cutoff and differ- 
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FIG. 9: Interaction measure A/T* in the QM model (upper panels) and PQM model (lower panels), see text for details. 



ent vacuum parameters. Interestingly, similar splitting 
effects also emerge in a finite volume FRG investigation 
within a two-flavour quark-meson model, cf . [72] . 

The splitting of the chiral transition is also reflected in 
a second minimum of the sigma meson mass as a function 
of the chemical potential as well as in the quark number 
susceptibility, which are sensitive to the chiral properties 
of the system. 

Increasing the pion mass towards its physical value we 
observe that the splitting in the chiral transition line per- 
sists. The outer branch of the transition is weakened, 
which can for example be seen in the behaviour of the 
sigma meson mass. At physical pion mass, the second 
transition is mostly washed-out. 

Of course, it would be interesting to study the vicin- 
ity of the critical point in more detail, but unfortunately 
some shortcomings of the Polyakov-loop potential pre- 
vent the computation of thermodynamic observables in 
the high chemical potential/low temperature region. The 
calculation of the phase structure itself relies mainly on 
the EoM, i.e., on derivatives of the effective potential 
with respect to the order parameters. These are not af- 
fected by the above-mentioned problems, which still al- 
lows us to study the phase structure itself. 

For a thorough study of the nature of the splitting re- 
gion at small pion mass, the deficiency of the Polyakov- 
loop potential poses serious problems. Regarding QCD, 
however, we are aware that in this high chemical poten- 
tial region our present approximation is incomplete. At 
baryon chemical potential [Ib = 3/i ^ 900 MeV, baryon 
degrees of freedom certainly play an important role in 
full QCD. Diquark fiuctuations and baryonic effects are 
expected to have a sizeable impact, as also suggested by 



two-colour results [64-67] and computations with iso-spin 
chemical potential, e.g. [73]. Thus, in order to produce 
reliable predictions of the physics in the vicinity of the 
QCD critical point - if it actually exists - an extension 
of the present model by the inclusion of these degrees of 
freedom is inevitable. 
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Appendix A: Differential Evolution Algorithm 

Usually, finding the roots of a higher-dimensional sys- 
tem of nonlinear equations, such as the EoM, Eq. (18), a 
natural choice would be Newton's method. This method 
converges quadratically if the initial values are already 
close to the final solution. Applied to our RG calcu- 
lation a huge number of iteration steps, of the order 
0(10'*), are actually needed to gain an acceptable nu- 
merical precision of the solution. However, each of these 
steps involves a complete RG evolution which increases 
the CPU time drastically. Hence, for the computation 
of thermodynamic quantities and higher derivatives, in 
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particular for low temperatures, an alternative method 
that requires less RG evolutions is inevitable. Such an 
alternative technique is provided by the differential evo- 
lution (DE) algorithm for global optimisation [74] which 
enables the solution of our three-dimensional system of 
the nonlinear EoM, Eq. (18), to high numerical precision 
with fewer RG evolutions. 

In detail the system Eq. (18) is solved as follows: we 
discretise the effective potential on a one dimensional 
grid (in radial direction a) and determine for arbitrary 
fixed initial values of the remaining variables ($, $) the 
potential minimum after the RG evolution. Then the 
EoM in a direction is fulfilled and we are left with a two- 
dimensional subsystem for the remaining Polyakov-loop 
variables. 

For these variables we define the cost function 

and apply the DE algorithm which we describe in the 
following. 

In the first step of the DE algorithm, an array Xij ofi = 
1, . . . ,iV pairs ($i, ^i), called target vector, is generated 
by randomly choosing values in a user-specified initial 



interval. For the present application we use ($i,$i) G 
[0, 1] X [0, 1] and typically the target vector has a length 
of TV = 10 - 20. 

Step two consists of defining the trial vector Vij , which 
is obtained from Xij by a randomisation procedure. For 
the details of this procedure, we refer the reader to [74]. 

Subsequently, the cost function is evaluated on the 
newly defined array Vij and compared to its value on 
the current target vector Xij. Note that this is the only 
step that involves the RG evolution. If, for a given index 
i, the result of the trial vector is smaller than the one of 
the target vector, the corresponding entry in the target 
vector is replaced by the one from the trial vector. The 
overall minimum of the cost function on the modified tar- 
get vector can then easily be calculated and compared to 
the user-specified criterion, e.g. /($min, ^"111111) < 10"'^. If 
this criterion is not fulfilled yet, the procedure is repeated 
from step two. 

With this algorithm it is possible to solve the EoM with 
an accuracy of at least 0(10"^) after O(IO^) generations, 
which corresponds to O(IO^) RG evolutions. For compar- 
ison, the same number of RG evolutions using Newton's 
method would typically yield roots of the order of 0(10^) 
only. 
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